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1 Abstract 

We investigate the dependence on length of op- 
tical fibres used in astronomy, especially the fo- 
cal ratio degradation (FRD) which places con- 
straints on the performance of fibre-fed spectro- 
graphs used for multiplexed spectroscopy. To 
this end we present a modified version of the 
FRD model proposed by Carrasco and Parry [1] 
to quantify the the number of scattering defects 
within an optical fibre using a single parameter. 
The model predicts many trends which are seen 
experimentally, for example, a decrease in FRD 
as core diameter increases, and also as wave- 
length increases. However the model also pre- 
dicts a strong dependence on FRD with length 
that is not seen experimentally. By adapting the 
single fibre model to include a second fibre, we 
can quantify the amount of FRD due to stress 
caused by the method of termination. By fitting 
the model to experimental data we find that pol- 
ishing the fibre causes more stress to be induced 
in the end of the fibre compared to a simple 
cleave technique. We estimate that the number 
of scattering defects caused by polishing is ap- 
proximately double that produced by cleaving. 
By placing limits on the end-effect, the model 
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can be used to estimate the residual-length de- 
pendence in very long fibres, such as those re- 
quired for Extremely Large Telescopes (ELTs), 
without having to carry out costly experiments. 
We also use our data to compare different meth- 
ods of fibre termination. 

2 Introduction 

Focal Ratio Degradation (FRD) is the non- 
conservation of etendue in an optical fibre, 
resulting in a broadening of the beam at the 
output. As noted by Parry [S], an increase in 
etendue within the instrument without a gain 
in information content is undesirable because 
the optical system which has to deal with it 
is consequently more complex and expensive. 
As telescopes increase in size, the large optical 
path lengths caused by the remoteness of focal 
stations from gravitationally-invariant instru- 
ment platforms and the increased sensitivity 
to mechanical flexure, mean that optical fibres 
will continue to be important for the key 
technique of highly-multiplexed spectroscopy 
|3j. Therefore, it is even more important to be 
able to predict the performance of optical fibres 
accurately in order to optimise the instrument 
system performance. 



In this paper we propose a modification to 
a model first proposed by Gloge [1] and later 
adapted by Gambling et al. [5] and Carrasco and 



Parry [lj. The model is described in section 3.1 



This modification aims to eliminate the depen- 
dence on FRD with length which is predicted 
theoretically but not seen experimentally. By 
modelling the fibre as two separate fibres with 
different amounts of scattering defects it can be 
shown that most of the FRD is caused by the 
stress frozen in at the fibre ends during end- 
preparation. 

The need for the fibre ends to be flat, smooth, 
and perpendicular to the axis of the fibre is 
well documented, although the method of end 
preparation varies widely depending on various 
restrictions. The quickest way to prepare a 
fibre is cleaving, typically by scoring it with a 
diamond pen or tungsten-carbide-edged blade, 
although a CO2 laser or spark erosion [7j 
has also been used. Once the fibre has been 
scored, it is placed under tension and then bent 
along the fracture until it snaps. Although 
this method is relatively quick to perform the 
results are not always satisfactory and a visual 
inspection must be performed. In addition to 
the problem of repeatability, it is extremely 
difficult to avoid a small defect occurring at 
the point where the fibre is scored as is shown 
in panel (a) of figure [TJ The second method is 
to polish the fibre on progressively finer grades 
of abrasive paper and finally on a solution of 
colloidal silica. This method of polishing the 
fibre is time-consuming, however visual and 
interferometric inspection of the end normally 
shows it to be high quality. Haynes et al. [8] 
have shown that the quality of the polish has 
a large effect on the FRD by comparing fibres 
with end face roughness of 245 nm and 6 nm rms. 



Figure [T] shows the typical finish which can 
be expected when a fibre is cleaved or polished. 
However all methods of preparing the ends of 
the fibre produce stress. The amount of stress 
and the depth to which it propagates from the 
end face is currently impossible to measure 
and methods can only be compared based on 
the smoothness of the finish as assessed with 
a microscope or interferometer. The most 
common rule-of-thumb is that stress propagates 
to 3 times the depth of the biggest defect on 
the end face of the fibre although this has never 
been proven. It has been shown by Craig et 
al. [9] that for the various fibres that they 
tested, there is no observable difference in FRD 
between a properly cleaved fibre and a polished 
fibre. 



In section 3.2 we modify the Gloge model to 
include two fibres with different properties in 
series. This allows us to model the end effect as 
a short fibre of fixed length with a larger number 
of scattering defects per unit length, do, than 
the main length of fibre. This parameter, do has 
different values depending on the termination 
method used. 

It has been shown by Avila et al. [10| that 
FRD is strongly affected by the way in which 
a fibre is mounted. Avila et al. tested UV 
grade silica fibres and found that, at an input 
focal ratio of F{ n = 8, F ou t fell from F/6.2 
when the fibre was held gently in a clamp to 
F/3.2 when it was squeezed strongly. This is 
an important issue because fibres are usually 
mounted in a ferrule or holder to protect the 
tip and allow them to be polished in bulk. 
When the fibres are mounted in this way they 
must be glued into position and it is extremely 
important to test the effects of various adhesives 
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Figure 1: Results of different methods of end- 
preparation for a typical fibre with a diameter 
of 100 \xm (a) unprepared (b) after cleaving (c) 
after polishing. In panel (b) the defect from scor- 
ing the fibre is clearly visible 



on the resultant FRD as shown by Oliveria et 
al. [TT]. This model also provides a quantitative 
measure of these mounting methods which will 
be independent of the fibre used. 



3 Theory Re-visited 
3.1 Gloge's model 

Carrasco and Parry [1] have shown that it is 
possible to measure the FRD of an optical 
fibre and characterise its performance using 
a single parameter, D. The model is based 
on a previous model by Gloge [I] who showed 
that the far-field distribution represents a 
direct image of the modal power distribution. 
Gloge developed a partial differential equation 
to describe the distribution of optical power 
P in a fibre of length L in terms of the ax- 
ial angle of incidence 6i n , with output angle 6 out . 



9 P (6 out, Gin) 

dL 



D__d 

din 9di 



A0 in P(0 O ut, 9 in) + 

dP(8 out An) 



(1) 



dOi: 



where A is an absorption coefficient and D is a 
parameter that depends on the constant cLq that 
characterises microbending: 



D 



A 



2dfn 



do, 



(2) 



A is the wavelength of light, df the fibre core 
diameter and n the index of refraction of the 
core. Equation [2] has been solved by Gambling, 
Payne, and Matsumure [5] for the case of a 
collimated input b e£tm of angle of incidence @in 
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P(0out,0in) = exp 

exp(-6L/2) 



1 — exp(— bL) 



Xi + X\ 1 + exp(-frL) 
1 — exp(— bL) 

(4 X g)5exp(-feL/2) 
1 — exp(— bL) 



(3) 



where X = {A/ D) 1/2 6 2 in , b = 4(AD) 1/2 , and 
Iq is the modified Bessel function of zeroth order. 

Equation [4] describes the angular dependence 
on the output flux in the far field for the case of 
a collimated input beam at an angle of incidence 
Bin . It is possible to generalise the solution given 
by equation [4] to the case of an input beam with 
any aperture as follows: 



F(6out, &0) 



Jo* Jo ' G(O in ,(f>,9o) 
x P{9out, 0i n )smdi n d9i r 



(4) 



where s'm6i n d9i n d(f> = dQ is the differential of 
solid angle, eft is the azimuthal angle measured 
around the same axis, and G is a function that 
represents the input beam. 

Carrasco and Parry have shown that when 
the value of D is found for a specific fibre, 
the results of other experiments where the far 
field output beam is recorded for input pencil 
beams of varying focal ratio, can be successfully 
predicted. By testing a 200 fim core fibre 
Carrasco and Parry found a solution to D which 
produced a reasonable fit to the experimental 
data. From their data, the accuracy of optical 
alignment can be inferred to be 1 deg since this 
shifts the predicted output focal ratio by 2.5% 
which is the typical size of their experimental 
errors. 



As shown by equation [2] the D parameter 
depends on the fibre core diameter, df, the 
microbending parameter, do, the wavelength 
of the injected light, A, and the length of 
the fibre, L. From this equation it is clear 
the D oc A 2 , and hence FRD will increase 
as the wavelength increases. This result has 
been verified experimentally [12J, although the 
opposite trend has been reported by Murphy et 
al. p2], who found a weak dependence on FRD 
with wavelength. However, measurements were 
made at F{ n = 3.5 where FRD is relatively 
unimportant compared to slower beams where 
the output beam approaches an asymptotic 
value. 

FRD be quantified conveniently by the FRD 
ratio Fi n /F ou t. The second trend predicted by 
equation pj is that D oc dj 1 . Figure |2 shows how 
this modeT and the two-fibre model aescribed in 
section [3T2] , predict how different parameters will 
affect the resulting output power distributions. 

Whilst the trends shown for the fibre core 
and do parameter dependence have been demon- 
strated experimentally as well as theoretically, 
the extent to which the length affects the FRD 
has not. Gambling et al. did report a length- 
dependence, however this result was for liquid- 
core fibres. This result has never been repro- 
duced for monolithic fused silica fibres. To con- 
firm this assumption an experiment was set up 
as described in section |4j The following section 
gives the theoretical description of the method 
of eliminating this length dependence. 

3.2 The two-fibre model 

In order to add a second fibre into the model 
it is assumed that for a single input ray at 
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output angle 



(a) (b) (c) 

df = 25 nm(solid) , 75 fj,m(dashed) , do = l(solid),20(dashed), L = l m(solid), 10 m(dashed), 

100 fj,m(dot — dashed) 50 (dot — da shed), 100 (dotted) 100m(dot — dashed) . 

Figure 2: Normalised theoretical results to show how the output power distribution will be affected by changing 
the input parameters. All results are for a solid input beam of angle 5° and unless otherwise stated for a 
fibre of length 1 m, df = 100 fim, and do = 1. 



angle 9i n to the fibre axis, the output surface 
brightness at angle 9 ou t is given by the Gloge 
model as F(9 out ,9 in ). 

For a solid input angle defined by unit surface 
brightness for 9i n < 6>o, the total output surface 
brightness is given by: 

Fi(9 ut,i, 9q) = / F(9 ut,i,9i n ^)2ir9i n< ic\9i n ^ 
Jo 

(5) 

If a second fibre is added to receive light from 
the first fibre, the equivalent number of single 
rays it sees at the input is given by: 

dN(9 out:1 ) oc Fitfout,!, 9o)2ire outtl de ou t,x, (6) 

so the output surface brightness in terms of 
the angle at the exit of the second fibre, # «t,2j 
is: 



F2(9 ut,2,9o) — / F(9 ou t,2, "out,! 
JO 

roc 

= / F(9 out ,2,9 ut,i)Fi(9out,i,9o) 
Jo 

x 2ir9 ou t,id9 out ,i 

(7) 

The coupling between fibres is achieved by 
setting 9in t 2 = 6<nit,i- 

Figure [3] shows curves for the output power 
distribution produced by both the one-fibre and 
two-fibre model. Since we can't determine do^ 
and L>2 independently, we use the product S = 
do, 2-^2 to characterise the second fibre. This pa- 
rameter essentially quantifies the total number 
of scattering centres (microbends ) present in the 
second fibre. The value of S was chosen to pro- 
duce the same amount of FRD at the asymp- 
totic value as the one-fibre model. This figure 
shows that the two-fibre model agrees well with 
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Figure 3: output power distributions predicted by the 
one fibre (dashed) and two fibre models (solid, with an 
offset of 0.03 in order to be able to show both curves). 



the predictions made by the one-fibre model (it- 
self rigorously tested against experimental data 
here and elsewhere). Indeed the difference in the 
95% enclosed energy values of the one- and two- 
fibre models differ by the same amount as the 
experimental error bars (see below). We use the 
angle which encloses 95% of the energy to define 
the output focal ratio since this reduces sensitiv- 
ity to stray light in the experimental setup. 



4 Experimental results 

To test the experimental dependence on FRD 
with length, a fibre was illuminated with a solid 
input beam of varying Bin as shown by figure 
[4j The fibre used was a typical fibre made by 
polymicro ( FIP100110125) with a core diameter 
of 100 fim, which has similar characteristics to 
fibres used in astronomical instruments. Input 
beams ranging between f/1.7 and f/18.9 (i.e. 
cone half-angles 16° — 1.5°) were used to in- 
vestigate how the FRD evolved with input angle. 

Tests were taken for a case of a cleaved fibre 
of length 10 m immersed via index-matching gel 
to a smooth flat glass plate and for the same 
fibre after it had been polished to a 1 fim opti- 
cal finish. The FRD curve presented in figure [5] 
show that the fibre has a fairly good FRD per- 
formance with an FRD ratio of around 1.1 at 
Fi n = 4.5. The rationale for testing a cleaved 
fibre with index matching gel was to make sure 
that we isolated the effect of stress from any scat- 
tering caused by imperfections in the topography 
of the end face. The use of index-matching gel 
for the cleaved fibre should eliminate the effect 
of surface roughness and the smoothness of the 
finish would do the same for the polished fibre. 
Furthermore, we expected that cleaving would 
result in relatively low stress and therefore pro- 
vide a baseline against which the results of the 
polishing could be compared. 

The test was repeated after the fibre had 
been cut to a length of 1 m for both the cleaved 
polished case. These two lengths of fibre will 
be referred to as 'long' and 'short' for the 
remainder of this report. 

The measured FRD for all scenarios under 
which the fibre was tested are also shown in fig- 
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Figure 4: Experimental setup. The adjustable iris was used to give a range of input focal ratios 
and a beam splitter was used in order to provide an image of the input light on the end face of the 
fibre. 



ure [5} From this, the D parameter for the fibre 
has been estimated in each of the specific cases 
(Table [I]). Theoretical predictions using these 
parameters are shown in figure [5] From this we 
estimate that the difference in FRD due to the 
different fibre lengths for the cleaved fibre at the 
asymptotic value was 0.15 ± 0.7, and hence is 
clear that within our experimental uncertainty, 
no length dependence is observed. 



5 Discussion of results 

We now use the experimental results to elimi- 
nate the length dependency and compare stress 
induced by different end preparation techniques. 
To this end the model can be run for various 
combinations of do,i and S (where S = ^0,2-^2) 
to find a single set of values which predict the 
experimentally-observed FRD curves. In order 
to find this set of values, a number of constraints 
can be applied to the two-fibre model, such that 
we require: 



1. do,i, S combinations where the difference in 
asymptotic F out for the long and short fibres 
is less than the experimental errors; 

2. do, i,5 combinations which give an F ou \, 
which agrees with the experimental data at 
the asymptote for both the cleaved and pol- 
ished fibres ; 

3. do 1 value less than that determined from 
the one-fibre model, and the same for both 
the cleaved and polished fibres. 

As is shown in section [4j the length depen- 
dence was less than i ? /0.15±0.67 for the cleaved 
fibre. The contour plot shown in panel (a) of 
figure [6] defines the upper limit on allowable 
AF, and all values within the 0.7 contour region 
give acceptable solutions for specific values of 
do,i and S . Secondly, the asymptotic -F out must 
be the same as the experimental data within 
the error bars. For the case of the cleaved fibre, 
the error in F ou i was determined to be 0.46 and 
this contour is shown in panel (b) of figure [6} 
All values within this 0.46 contour region give 
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Table 1: Parameters need to produce theoretical curves which match experimental results for the 
fibre preparation methods. 



preparation 


fibre 


FRD ratio 


do 


D 


method 


length 


(at F in = 4.5) 


(m _1 ) 




cleave 


long 


1.07 


5.0 


2.1 x 10" 5 


cleave 


short 


1.09 


52 


2 x 10~ 4 


polish 


long 


1.19 


9.1 


4 x 10~ 5 


polish 


short 


1.15 


85 


3.7 x 10~ 4 




Figure 5: FRD curves where experimental results are shown by crosses and theoretical results 
produced by the original one-fibre model are given by a smooth line 
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acceptable solutions for specific c?o,i and S . 
The final panel in figure [6] shows how these con- 
straints combine, and a solution for the cleaved 
fibre is shown by a cross at do,i = 0.1 and 
S = 25. Using these constraints for the polished 
fibre produces a solution at do,i = 0.1 m _1 and 
5 = 48. 

Figure |7] shows the FRD curves for various 
lengths of fibre using these parameters. From 
this figure it is clear that we have removed 
the length dependence of the model within 
the lengths for which it was constrained. The 
small discrepancies between the theoretical 
and experimental results arise due to mis- 
alignments within the experiment and the 
fact that we have assumed G was a perfect 
step function rather than measuring it explicitly. 

The second objective of the two fibre model is 
to provide a method of comparing stress induced 
by various end preparation techniques. Compar- 
ison of the parameter, S, between the polished 
and cleaved cases gives a quantitative assessment 
of the induced stress. The relative values show 
that there are twice as many microbends per 
unit length when polishing the fibre, compared 
to when it is cleaved. The high absolute values of 
S make a very strong argument for taking great 
care to eliminate stress in the end-preparation 
process. 

Currently, there is no method of determining 
the depth to which stress in the end of the fi- 
bre extends. The rule-of-thumb that the stress 
propagates to 3 times the depth of the largest de- 
fect on the end face of the fibre implies that the 
stress only affects a few //m of the length of the fi- 
bre and forces do, 2 to be extremely large. Thus, 
assuming that the defect size is determined by 
the smallest grit size, implies L2 = 3 \im and 



do,2 = 8.3 x 10 6 mT 1 and d ,2 = 16.0 x 10 6 m _1 
for the polished and cleaved fibres respectively. 

6 Conclusions 

We have adapted the power distribution model 
described by Gloge and produced a new model 
which will eliminate the length dependence on 
FRD in order to provide a better agreement 
with experimental data. This model has been 
tested against the original one-fibre model, and 
with experimental data, and shown to be in good 
agreement. As a result of trying to quantify the 
amount of stress within the end of the fibre it 
has been shown that the end-effect is extremely 
powerful. If we assume that the do,2 parameter 
characterizes the FRD induced in the end of the 
fibre then it is possible to compare how the FRD 
produced in the majority of the length of the fi- 
bre is affected. In the case of the cleaved fibre 
the do parameter (which quantifies the amount 
of microbending) calculated using the one fibre 
model was reduced from 5 m^ 1 for aim fibre 
and 52 mT 1 for a 10 m fibre to 0.1 m^ 1 indepen- 
dent of length using the two-fibre model. This 
significant reduction is testament to the strength 
of the end effect and clearly shows how much 
FRD is generated in the end of the fibre. 

It was initially assumed that cleaving the fibre 
would produce very little stress, however this was 
not shown to be the case. This model is there- 
fore extremely useful when comparing various 
end preparation techniques. For example, our 
results show that cleaving the fibre halves the 
strength of the end effect compared with polish- 
ing. The model can also provide an upper limit 
on the FRD performance of any length of fibre 
and will therefore prove to be very useful when 
building instruments for Extremely Large Tele- 
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(a) difference in asymptotic F ou t between L=l m 
(b) difference in asymptotic F ou t between experimental 
and L=10 m with contours shown in increments of 0.7 and theoretical results for the cleaved fibre 
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Figure 6: The two fibre model was run for various combinations of do,i an d 5* (where S = do, 2-^2) • 
Panel (a) shows the difference in asymptotic F ou t between L=l m and L=10 m with contours shown 
in increments of 0.7. Panel (b) shows the difference in asymptotic -F ou t between experimental and 
theoretical results for the cleaved fibre. Panel (c) shows how these graphs can be combined to find 
a single solution, marked by the large cross, for c?o,i and S which satisfies experimental data. The 
black vertical line in panel (c) shows the do,i parameter predicted by the one fibre model for the 
cleaved fibre. 
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Theoretical FRD curves 
2 fibre model - cleaved fibre 



Theoretical FRD curves 
2 fibre model - polished fibre 




Figure 7: 2 fib model results with do,i = 0.1 m and ^0,2-^2 = 25 for the cleaved fibre and 
^0,2-^2 = 48 for the polished fibre. It is clear that within the errors the length dependence of the 
original model has been removed. 



scopes. 
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